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1 INTRODUCTION 



ABSTRACT 

The alignment of dark matter (DM) halos and the surrounding large scale structure 
(LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital 
angular momentum of subhaloes is investigated relative to the LSS using the eigen- 
vectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep 
within the non-linear regime. Knots, filaments, sheets and voids are associated with 
regions that are collapsing along 3, 2, 1 or principal directions simultaneously. Each 
halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to 
the nature of the collapse at the halos position. The full distribution of shear eigenval- 
ues is found to be substantially different from that tagged to haloes, indicating that 
the observed velocity shear is significantly biased. We find that larger mass haloes live 
in regions where the shear is more isotropic, namely the expansion or collapse is more 
spherical. A correlation is found between the halos shape and the eigenvectors of the 
shear tensor, with the longest (shortest) axis of the halos shape being aligned with the 
slowest (fastest) collapsing eigenvector. This correlation is web independent, suggest- 
ing that the velocity shear is a fundamental tracer of the halo alignment. A similar 
result is found for the alignment of halo spin with the cosmic web. It has been shown 
that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at 
which this spin flip occurs is web dependent and not universal as suggested previously. 
Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, 
providing a possible insight into the highly correlated co-rotation of the Milky Ways 
satellite system. The present study suggests that the velocity shear tensor constitutes 
the natural framework for studying the directional properties of the non-linear LSS 
and of halos and galaxies. 



attempt to quantify the growth of perturbations and thus 



According to the current cosmological paradigm, the early 
Universe was seeded with small density perturbations that 
peppered an otherwise homogeneous matter distribution. As 
the Universe expanded, these peaks in the density field grew 
due to gravitational instability at the expense of under dense 
regions producing a highly complex network of voids, walls, 
filaments and knots, known as the cosmic web. The first 



the cosmic web is attributed to |Zel'dovich ( 1970 1 who's fa- 
mous "pancake" solution forms the basis of how we think 
about the cosmic web. Further work (by, e.g 



Bond et al. 



[19961 [Shandari n et al.|2004l|Shen et al.|2006[ and references 
therein) began to examine in great detail the nature of how 
the cosmic web formed and its defining properties. 

This theoretical view is corroborated by mapping the 
distribution of galaxies on the largest scales (by e.g. the 
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Figure 1. The cosmic web in a 0.97 Mpc/h (one grid cell) slice. Left panel: the DM density distribution, with contours denoting regions 
that are over dense with respect to the mean. Middle panel: The cosmic web in the same slice with the most massive haloes over plotted 
in red. The hierarchy of cosmic web types is visible as knots (black) are embedded in filaments (light gray) which are embedded in sheets 
(dark gray) which surround voids (white). Note that due to the thinness of this slice, sheets appear one dimensional. We over plot regions 
above the mean over density as empty contours. Right panel: A zoom into a 50 Mpc//i sub region (denoted in the middle panel by the 
small square) showing each halo as an ellipse whose major and minor axes are defined by the inertia tensor. Again, regions more over 
dense than the mean are shown with the thin black contours. The haloes are aligned with large scale structure as shown in Fig. [6] Note 
that the size of each halo is not representative of the halo's virial radius. 



2dFGRS and SDSS, |Colless et ai^|200l) |York et al.|[2000]) 



which quite clearly defines the cosmic web (Erdogdu et al 



20041. Many studies have examined web attributes in de- 
tail fo cusing on knots (jEinasto et al.|1997|), filaments ( Coles| 



1996 Bharadwaj et al. 20041, sheets (Einasto et al. 20111 



and voids ( Kreckel et al.|20"IT Pan et al.|2012 Cuesta et al 
20081. Recently, Munoz-Cuartas et al. (20111 used the data 



from the SDSS to reconstruct the cosmic density field, in- 
cluding the cosmic web. 

Identifying the cosmic web in numerical simulations 
poses a number of interesting challenges. The simplest tech- 
niques (e.g Avila- Reese et al.|2005 Sousbie et al.||2008 \ de- 
fine the web purely geometrically, often only in terms of the 
density field smoothed over a few Megaparsecs. |Hahn et al.| 
( 2007 1 and Forero- Romero et al. ( 2009 1 developed the con- 



cept using a dynamical stability criterion (the Hessian of the 
potential) which in many ways is more fundamental than the 
density field since it reflects the dynamics of matter. Since 



the tidal web of Hahn et al. ( 2007 1 is in principle only ap- 



plicable to regions where the tidal tensor is representative 
of the large scale dynamics, it is incapable of identifying 
the finer web seen in high resolution numerical simulations. 
Hoffman et al. (20121 extended these ideas by suggesting 



that the velocity shear (so-called "V-web") is a better rep- 
resentation of the cosmic web on linear and sub-linear scales 
(the velocity shear tensor converges to the tidal tensor in the 
linear regime). 

While the cosmic web is interesting in and of itself, its 
characteristics affect the haloes and galaxies that inhabit 
it. For example, tidal torque theory (Peebles 1969 White 



1 1984[ > , the theory by which galaxies acquire angular momen- 
tum, suggests that the large scale tidal force field is respon- 
sible for both the formation of the cosmic web and the spin 
of galaxies. A correlation between galaxy spin and web en- 
vironment is thus expected. This correlation has been found 



in a few observational studies (Binggcli 1982 Tempel et al 



2012 1 that defined aspects of the cosmic web using various 



techniques. The recent work of Tempel et al. ( 2012 1 found a 
weak (but significant) trend for the spin axis of disc galaxies 
to be aligned with filaments, while that of ellipticals appears 
to be perpendicular to them. Yet these observational stud- 
ies of environment are really in their infancy, since obtaining 
a three dimensional spin vector for each galaxy as well as 
a proper classification of the cosmic web is difficult (if not 
impossible when, for example, considering the angular mo- 
mentum of elliptical galaxies). 



A number of numerical studies (e.g. 


Altay et al.||2006 


Aragon-Calvo et al.||2007| |Brunino et al.| 


2007| |Codis et al.| 


2012 


) have examined the orientation of halo shape or angu- 



lar momentum with respect to the cosmic web in just one or 



two given web types. For example Altay et al. ( 2006 1 looked 
at the alignment between halo shape and the main filament 
axis while Codis et al. (20121 examined the orientation of 



halo spin for haloes found in filaments. | Aragon-Calvo et al.| 
(20071 looked at the orientation of angular momentum of 
filament and sheet haloes while Brunino et al.l (120071) ex- 



amined both spin and shape alignments in voids found in 
the Millenium simulation ( Springel et al.|[2006 1 . The recent 
work of Trowland et al. ( 2012 1 also used the Millenium sim- 
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Figure 2. The probability density distribution for the largest (left), intermediate (middle) and smallest (right) eigenvalues of the velocity 
shear tensor. The distribution of eigenvalues for all cells is shown in green. We also present the distributions of eigenvalues assigned to 
low (black), intermediate (blue) and high (red) mass haloes. Each distribution is normalized to the number of objects it contains. The 
noise in the high mass (red) curve, is due to poorer number statistics in that mass bin. Note that the distribution of eigenvalues assigned 
to haloes are significantly different from those of the full cosmological volume: just as haloes trace the highest peaks in the density field, 
they are also a biased tracer of velocity shear. 



ulation to examine the evolution of low and high mass halo 
spin alignments within filaments. 

The general consensus from most of these studies is two 
fold: the first is that halo shape correlations are stronger 
than spin correlations. Second is that haloes are oriented 
(at the ~ 20% level) with the cosmic web. Haloes tend to 
point along the spines of filaments, parallel to the plane of 
cosmic sheets, and perpendicular to the radial direction of 



voids. |Arag6n-Calvo et al.| ( |2007| ) and |Codis et aL] ( j2012 | 
found that the angular momentum axis of low mass haloes 
points along the filament axis while that of high mass haloes 
points perpendicular to it. The mass at which this spin flip 
is seen was reported to be ~ 1 — 5 x 1O 12 M0 and is possi- 
bly due to the fact that high mass haloes grow by mergers 
which preferentially come from the filament axis, while low 
mass haloes are "wrapped up" with the filament as it col- 
lapses. While clearly an attractive explanation, the critical 
mass of the alignment transition appears empirical and ar- 
bitrary since it has not been theoretically motivated. More- 
over, Trowland et al. (20121 examined the evolution of the 
spin alignment with redshift and found a completely op- 
posite effect: all haloes are oriented perpendicular to their 
filaments axis at high redshift and it is only low mass haloes 
which become parallel later on. Clearly the origin of the spin 
web alignment requires significant study. 



Libeskind et al. (20121 used the V-web to examine the 
orientation of haloes and, for the first time, their subhaloes 
with respect to the cosmic web in simulations tuned to re- 
produce the local universe. They found that the V-web con- 
firmed previous studies that showed that high (low) mass 
halos tend to spin perpendicular (parallel) to the defining 



filamentary axis. Due to its limited sample size it was un- 
able to measure the transitional mass at which halo spin's 
flip. However, they found that subhalo orbits as well as halo 
spins are correlated with the large scale velocity shear. This 
result was surprising since subhaloes orbit deep in their par- 
ent potential, in a way shielded from the large scale velocity 
shear. That the cosmic web can still be reflected amidst 
subhalo orbits, is a testament to its importance in galaxy 
formation. 

The finding that subhaloes reflect the cosmic web is 
suggested by studies that focus on the alignment of subhalo 
positions and velocities with respect to each other and their 



host (|Libeskind et al. 2005 


Zentner et al.||2005 Vera-Ciro 


et al.||2011 Libeskind et al. 


2011 e.g.). Yet Kroupa et al. 



gued that the satellites of the Milky Way can not be ACDM 
substructures since their anisotropic spatial and kinematic 
z = distribution is incompatible with simulations, which in 
a naive sense predict an isotropic spatial distribution with 
randomized orbits. Thus the alignment between substruc- 
tures and their large scale structure is of great interest in 
the study of the Milky Way. 

In this paper, we extend the work of Libeskind et al.] 



(2012) and examine halo and satellite alignments with re- 
spect to the large scale structure as defined by the V-web. 
We take a detailed look at both the cosmic web itself, by 
studying the eigensystems of the diagonlized shear tensor, 
and the effect on haloes within it. Numerical methods and 
the algorithm used to dissect the cosmic web are explained in 
Section [2] The main results with respect to the cosmic web 
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and halo alignment are presented in Section [3] We conclude 
in Section. [4] 



Fractional Anisotropy 



2 METHODS 



2.1 The Bolshoi Simulation 



In order to investigate correlations with the cosmic web, 
we use the Bolshoi ( Klypin et al.|20TT l Af-body simulation. 
We refer the reader to |Klypin et al.| ( |2011[ ) for details on 
the simulation, and include just the salient points here. The 
run follows the evolution of matter in a large periodic box 
according the ACDM paradigm, employing cosmological pa- 



rameters that are compatible with WMAP7 (Jarosik et al 



2011[ ) data, e.g. h = 0.7, ST2 A = 0.73, fi m = 0.27, Q b = 0.046, 
n = 0.95 and as = 0.82. Dark matter (DM) is numerically 
simulated using an Adaptive-Refinement-Tree code (ART, 
Kravtsov et al.|1997 |. The simulation box has a side length 
of 250 ft" 1 Mpc ■ and is filled with 2048 3 w 8 billion particles, 
each with a mass of ~ 1.3 x 10 8 /i _1 M©. The proper physi- 
cal spatial resolution (smallest cell size) of the simulation is 
lft -1 kpc. In this paper we use only the simulations z = 
snapshot. 

DM halos are identified using a Friends-of-Friends (FOF 
Davis et al.|1985| ) algorithm with b = 0.17. FOF links parti- 
cles closer than b times the mean inter particle separation, 
into objects referred to as haloes. In effect, the free parame- 
ter b defines isodensity contours in the matter distribution. 
The choice of b — 0.17 returns haloes whose mean over- 
density is roughly equal to that of virialized objects in the 
CDM cosmology. For each halo, internal properties such as 
mass and angular momentum (spin), are calculated. The 
shape of each halo is defined by diagonalizing the reduced 
| |Gerhard|1983 Bailin fc Steinmetz|2005 \ inertia tensor hj , 
constructed by summing over n particles in the halo: 



E 



(i) 



where r n = (x^+y^+z^) is the distance to each nth particle. 
The reduced inertia tensor is used instead of the standard 
one (namely without the r\ term in the denominator) so as 
to lessen the effect of recently accreted large subhaloes in the 
outer parts of the halo. The three eigenvalues (a, b and c) are 
sorted in increasing order such that a > b > c. The halo's 
principle axes (a, b and c) are associated to the eigenvectors 
that correspond to each eigenvalue. An arbitrary mass cut 
of 1O 9 ' 5 /i -1 M0 is applied to the halo catalogue such that 
very small haloes are excluded from the analysis. 

Substructures are identified using the Bound-Density- 
Maxima BDM technique (Klypin fc Holtzman|1997[ ). BDM 
finds local maxima in the density field, removes unbound 
particles, and computes the relevant halo properties (such 
as mass) within a truncated radius. We impose a 30 parti- 
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Figure 3. The Fractional Anisotropy (FA), as defined in equa- 
tion [6] is a measure of shear anisotropy. The same slice of the 
simulation shown in Fig. [T] The thin white lines separate voids 
from sheets, filaments and knots. The continuous black and white 
regions denote low and high values of FA, ranging from zero to 
unity. Note how the higher values of FA are found abutting the 
contours and snaking through the voids, revealing the inner struc- 
ture of voids. 



cle limit on the BDM subhaloes, ignoring smaller clump^] 
Note that the BDM catalogues are only used to identify the 
orbital angular momentum of subhaloes. 

Throughout this paper, haloes are grouped according to 
mass into three bins, a low mass bin (M V - 1T < 10 11,5 /i _1 Mq, 
colored black), an intermediate "Milky Way" mass bin 
(1O 11 ' 5 /i" 1 M < Mvir < lO 12 - 5 ft _1 Af , colored blue), and 
a high mass bin (Af vir > W 12 - 5 h~ 1 M e , colored red). 



2.2 The Cosmic Web 

In order to quantify the cosmic web, the velocity shear ten- 



sor of the particle distribution is computed (see |Hoffman| 



et al 



2012 for a detailed analysis of this technique). First, 



a 256 grid based density and velocity field is constructed 
using the "clouds in cells" (CIC) technique, resulting in a 

1 The Bolshoi simulation (halos as well as the DM particle po- 
sitions and velocities) can be downloaded from the MultiDark 
database (http://www.multidark.org). The databases as well as 
the FOF and BDM halo finder used here are described in IRiebel 
|et al.|(2uTf) . 



© 2002 RAS, MNRAS 000, ??-?? 



The velocity shear tensor 5 



< 

o_ 
>-, 



O 




Finding the threshold which best returns the visual im- 
pression of a cosmic web is not self-evident, as it depends 



0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 4. The Fractional Anisotropy FA is a measure of shear 
anisotropy. Since each halo can be assigned a web type based 
on the local velocity shear calculated at its location, we divide 
our halo sample into four web classes, plotting the probability 
density distribution of FA for haloes resident in knots (black), 
filaments (green), sheets (blue) and voids (red). Each distribution 
is normalized by the number of objects in it. FA=0 indicates 
uniform expansion or collapse. Since the peaks and widths of these 
distributions depend on web type, we infer that each web type 
behaves significantly different: knot collapse is not simply void 
expansion played in reverse. 



smoothed density and velocity distribution. The CIC'd ve- 
locity field is then Fourier Transformed into fc-space and 
smoothed with a Gaussian of kernel width equal to one grid 
cell (0.97ft _1 Mpc) in order to wash out artificial effects 
introduced by the preferential axes of the Cartesian grid. 
Using the Fourier Transform of the velocity field the nor- 
malized shear tensor is calculated as: 



1 

~2H~o 



dv a 



dvp 



drp dr a 



(2) 



where a, j3 = x, y, z. Similar to the inertia tensor, the eigen- 
values (Ai, A2 and A3) and corresponding eigenvectors (ei, 
e.2 and £3) of the shear tensor are obtained at each grid cell. 
Following convention, the eigenvalues are ordered such that 
Ai > A2 > A3 and a web classification scheme based on how 
many eigenvalues are above an arbitrary threshold is car- 
ried out. If none, one, two or three eigenvalues are above 
this threshold, the grid cell is classified as belonging to a 
void, sheet, filament or knot. 



on the numerical parameters used. |Hoffman et al.| ( 2012 1 
and Libeskind et al.N2012|), found a threshold of A t h =0.44 



best reproduced the visual impression of the cosmic web. 
This choice of threshold implies that when the expansion 
(or contraction) is very weak, it is considered contraction 
(or expansion). This threshold is used in this work as well, 
following the same visual justification. We defer to a future 
study the effect of varying the threshold on the cosmic web. 

The cosmic web of the simulation is presented in Fig [I] 
where we show a thin slice through the computational box. 
The density field is depicted in the left panel, with contours 
enclosing regions that are above the mean density. This can 
be compared with the cosmic web of the same slice shown in 
the middle panel of Fig, [I] (with haloes over plotted in red). 
The reader will note that more than 80% of the regions of 
the simulation that are above the mean over density are des- 
ignated as sheets, filaments and knots by the web classifier. 
The right most panel of Fig, [T] shows a zoom on a random 
sub-volume in the box, chosen to highlight the alignment 
of haloes with the large scale structure, as described in sec- 
tion E21 

Haloes (and their subhaloes) are connected to a local 
velocity shear tensor by simply identifying the CIC grid cell 
in which the halo sits. Properties of the cosmic web (as cal- 
culated by the local shear tensor) in that cell are then as- 
sociated to each halo within it - specifically, each halo is 
assigned environmental attributes, such as shear eigenval- 
ues, eigenvectors and web type. In Section [3. 1| we associate 
the uniformity of the velocity shear to haloes as well. 

The discrete sampling of the cosmic density field is an 
attempt to numerically represent the underlying fluid. Natu- 
rally, any technique which attempts to redraw a continuous 
field in discrete terms will undoubtedly be plagued by ar- 
tifacts introduced by the sampling itself. When examining 
the cosmic density field, these problems are acute since the 
density field spans many orders of magnitude and is consid- 
erably heterogenous, while the grid is regular. Ideally the 
cell size and smoothing kernel used would be adaptive and 
varied based on the local properties of the underlying den- 
sity field. In the absence of an adaptive grid, the CIC and 
the Gaussian kernel chosen for the smoothing set a physi- 
cal scale on the objects which can be examined. A grid size 
that is smaller than a halo's virial radius will probe the shear 
tensor interior to that halo, whereas a much coarser one will 
calculate the shear in a larger region, thereby washing out 
any inherent spin- web correlation. Thus the careful selection 
of an appropriate grid size and smoothing scale are crucial. 

In this work we choose a 256 3 grid, which spans the 
250Mpc/ft cosmological box, giving a cell size of « lMpc/ft. 
All objects whose virial radius is smaller than this will be 
embedded in a cell: some will fit snugly, others more loosely. 
This grid size has been chosen such that all the objects in 
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our halo catalogue can be reliably examined (i.e. a negligi- 
ble number of haloes have r v ir > lMpc//i.) In practice this 
choice of grid puts an upper limit of ~ 2 x 1O 14 M0 jh on the 
halo mass we may probe. 



3 RESULTS 

Before presenting results on the alignments of haloes with 
respect to the eigenvectors of the velocity shear tensor, we fo- 
cus on understanding the distribution of the eigenvalues, Ai, 
A2 and A3 . The range of the smallest and largest eigenvalues 
is essentially bounded by the value of A2 , while the interme- 
diate eigenvalue is constrained by the values of Ai and A3. 
That is, 

Ai G [A 2 ,oo) 
A 2 G [A 3 ,Ai] 
A3 G (-00, A 2 ] 

The eigenvalue describes the strength of the expansion 
or collapse of matter along the corresponding eigenvector. 
For example, three large positive eigenvalues indicate strong 
collapse along all three axes. On the other hand in filaments 
and sheets where Ai> A t h and A3< A t h expansion along 
one axis is roughly equal to collapse along another when 
A3 I ~ |Ai|. Note that the trace of the shear tensor is 



Tr(E) = 



Ai + A 2 + A3 
= -V-Voc<5 



(3) 
(4) 



and is thus proportional to the over density 8 in the linear 
regime. 

We now turn to the distribution of velocity shear eigen- 
values and show, in Fig. [2j the normalized probability den- 
sity distribution of Ai, A 2 , and A3. Since each halo can be 
assigned an eigen-system based on its location, both the full 
distributions of eigenvalues for all cells (green lines), as well 
as subsets of eigenvalues that are assigned to haloes binned 
by halo mass (black, blue and red lines) are shown. Note 
that the latter is identical to the distribution of eigenvalues 
weighted by halo number. Each distribution is normalized 
to unity. 

The full distribution of A's (green curves) shows a num- 
ber of interesting properties. Firstly, they are all well repre- 
sented by log-normal distributions, whose median and width 
are eigenvalue dependent. Secondly, A3 is not simply a neg- 
ative reflection of Ai, as would be expected if expansion 
and collapse in different regions of space directly correlated 
with each other (i.e. if collapse looked like expansion in re- 
verse). Instead the width of the distributions decrease with 
decreasing eigenvalue: Ai has a relatively wide width while 
A3 is much narrower. This implies that expansion or collapse 
along ei may occur with varying strengths while the range 
of the expansion or collapse strength along £3 is more tightly 
restricted. 



The distribution of eigenvalues for those regions of space 
that contain haloes (black, blue and red curves in Fig. [2| 
shows a number of interesting attributes as well. Despite 
all being drawn from the same parent distribution (green 
curves), they are qualitatively very different. Firstly we note 
that none of the eigenvalue distributions (when binned by 
mass) resemble the parent distribution in terms of width or 
median. The low mass bin (black curve) peaks at a similar 
value to the parent distribution, probably owing to the fact 
that this mass bin spans the largest fraction of grid cells. Yet 
the peak of each distribution is clearly mass dependent. The 
more massive the set of haloes being examined, the greater 
the median of the distribution. Intuitively this is because the 
most massive haloes reside in regions of strong gravitational 
collapse that is, regions with larger values of A. 



3.1 The Fractional Anisotropy: The distribution 
of shear sphericity 

The semi-infinite range of eigenvalues (described above), 
makes it difficult to characterize the uniformity of the col- 
lapse or expansion in the way that, for example, the inertia 
tensor may be characterized by its triaxiality. In their sem- 



inal work on Gaussian random fields, Bardeen et al. ( 1986 1 



defined the "eccentricity" e of a density field (not a shear), 



Ai - A 3 
2\6\ 



(5) 



where in the linear regime the dimensionless over density is 
defined as the trace of the shear, i.e. equation [3] This mea- 
sure of eccentricity has been used by a number of authors 
(e.g. |Lee fc Lee|2008l ) when examining velocity shears. That 
said, the original Bardeen et al. definition was invoked in a 
system with three positive eigenvalues (such as the inertia 
tensor) and not in the situation where the eigenvalues are 
bound by (— 00, c») as in the shear tensor. 

The usage of the Bardeen et al. eccentricity is there- 
fore highly problematic in our situation. In cases where 
Ai ~ —A3, small values of A 2 can cause the denominator 
of equation [5] to tend to making e highly unstable. Also, 
the fact that e is semi-infinite makes it more difficult to 
compare vastly different eccentricities. 

In order to overcome these difficulties, the amount of 
anisotropy in the shear tensor can be characterized using 
the fractional anisotropy, FA 



FA = 



(Ai - A3) 2 + (A 2 - A3) 2 + (Ai - A 2 ) 2 



Aj + A 2 + A 3 



(6) 



which measures the fraction of the magnitude of the shear 
that is due to the anisotropy of the eigenvalues ( Basser 1995 1 
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[^] Note that the l/y/3 is introduced to normalize the frac- 
tion to unity. This definition is employed to characterize the 
uniformity of the collapse or expansion. 

The FA measures the kinematical "morphology" of the 
velocity shear. It is superior to the eccentricity since it is 
forced to be between zero and unity. FA=0 implies isotropic 
expansion or collapse in all three directions. In filaments 
and sheet volumes, FA=1 indicates that collapse along ei is 
roughly the same strength as expansion along £3, while in 
knots and voids, FA—1 indicates large anisotropics in the 
expansion or contraction. 

In Fig. [4] we show the distribution of the local velocity 
shear FA. Since each web type represents a different number 
of positive eigenvalues we plot them seperately. A number 
of important points can be gleaned from this figure. 

• The collapse in knot environments is relatively 
isotropic; that is, the value of FA measured in knots tends 
to be low. The strength of the velocity shear in knot envi- 
ronments thus displays a remarkable regularity. 

• The distribution of FA in filaments and sheets is wide - 
the shear that defines these web types takes on a myriad of 
kinematical morphologies. Whereas filaments peak at lower 
values of FA, sheets peak at higher values., 

• Void regions have high values of FA indicating highly 
anisotropic expansion. 

We thus conclude from Fig. [4] that the cosmic web re- 
flects a remarkable variety of shear anisotropies. Voids for 
example, are often not kinematically isotropic objects. Knots 
display significantly more uniformity than voids, although 
also span a wide range of sphericities. Filaments and sheets 
on the other hand, while displaying some similarities are still 
coloured by a variety of corpulence: some filaments are fat, 
some are thin, some sheets are thick some are narrow. 

The FA is shown graphically in Fig. [3] Since voids dom- 
inate the volume, the figure emphasizes the rich diversity in 
the kinematic "morphology" of voids. Although it appears 
that much of the void volume is represented by regions of 
high FA (white areas in Fig.|3|, void regions that abut sheets 
have lower FA (black regions in Fig. |3J. Interestingly this 
plot shows exactly why the FA spans the full range of values 
in voids: deep inside voids, structures of low FA can be seen 
snaking through the expanse. The FA is thus a measure that 
can be used to probe the inner morphological structure of 



voids similar to Aragon-Calvo & Szalay (2012|. Recall that 
the FA is insensitive to the sense of the anisotropy and will 



2 The concept of the Fractional Anisotropy was developed in 
the field of Nuclear Magnetic Resonance by |Ba sscr (1995J as a 
tool in brain imaging. The diffusivity of water molecules through 
cerebral tissue, one of the quantities measured by an NMR scan, 
constitutes an identical process to the velocity shear. Thus both 
the cosmic velocity shear tensor and the diffusion tensor in brain 
scans can be quantified by the Fractional Anisotropy. 




Figure 5. Schematic illustration of the angles examined in this 
work. The (reduced) inertia tensor is characterized by three prin- 
ciple axes (a, b, and c) that represent the long, intermediate and 
short axes of the halo respectively. The velocity shear is tensor is 
defined by the orthonormal basis ei, &2 and (=3, where the eigen- 
values and eigenvectors have been ordered to represent the axes 
along which material is collapsing fastest to slowest. We only con- 
sider the angles formed between vectors of the same line style, in 
other words c ■ e\ , b ■ ei and a ■ ez 



return the same value for uniform expansion as for uniform 
contraction. 

Armed with an understanding of the distribution of 
shear tensor anisotropy, we proceed to examine how haloes 
and their subhaloes are aligned with the cosmic web and 
how the kinematical nature of their web environment affects 
these alignments. 

3.2 Halo alignments with respect to the cosmic 
web 

The alignment between haloes and the large scale structure 
can be quantified by examining the distribution of angles 
formed between each eigenvector of the diagonalized shear 
tensor and a vector defining the halo. In what follows we 
investigate the alignment of the halo spin axis (J) and the 
halo shape (a, b, c) with the local velocity shear. Since FOF 
haloes are not rotationally, these alignments are not identi- 
cal. 

When examining the alignment of haloes and subhaloes 
with the principle axes of the velocity shear tensor we choose 
to divide the halo sample by web element and by mass. The 
number of haloes in each mass bin and web type is presented 
in Tabled] 
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Figure 6. The probability distribution P(|cos0|) as a function 
of the angle formed, |cos</>|, between the shape of a halo and 
the cosmic web. We show cos0 =e^-a (left column), cost/) =ei-b 
(middle column), and cos <f> =e\-c (right column). The full angular 
distribution function (top row) is complemented by examining 
the alignment as a function of web environment, that is in haloes 
found in knots (second row), filaments (thtird row), sheets (fourth 
row) and voids (bottom row). Haloes arc binned by mass, with 
low mass haloes (M < 10 11 h~ ^Mq) in black, intermediate mass 
haloes (10 11 - 5 h- 1 M Q < M < 10 12 ' 5 /i-lAf o ) in blue and the 
most massive haloes (M > 10 12,5 h~ 1 Mq) in red. The strength of 
each distribution is characterized by the average offset between 
it and a random distribution in units of the Poisson error and is 
indicated by the colored numbers at the bottom of each panel. 
Note that even though the low mass bins often look "flatter" , 
the signal is often more significant since the larger sample size 
decreases the Poisson error. The error bars in the left hand column 
indicate the expected poisson noise for a uniform distribution of 
that size. 
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Figure 7. The median (cosine of the) angle formed between 
£3 and the shape (top) and the spin (bottom) of haloes as a 
function of web type. The median angle for high, intermediate and 
low mass haloes are shown in red, blue and black respectively. A 
uniform distribution has a median of < cos</> = 0.5 >. 



3.2.1 Alignments with respect to halo shape 

In practice since there are two ortho-normal bases (the ve- 
locity shear tensor and the inertia tensor) there are a total 
of 9 (=3 2 ) possible angles to examine. Since many of these 
are complimentary we choose to examine the angles formed 
between ei and c, e.2 and b, and £3 and a. These angles are 
chosen since if, e.g. a prolate halo points along the filament 
it lives in, cos(f> = e3.fi, = 0. Similarly if a disc shaped halo 
lies in a sheet cos cj) = ei.c = 0. Note that since the principle 
axes of both tensors have no direction (that is, they define 
an axis, not a direction), only the absolute value of the angle 
is considered. 

In Fig. [6] we show the alignment of halo shape with the 
cosmic web. In the upper row we show the alignment for the 
full population, while in the second, third, fourth and fifth 
row, we show the alignments in knots, filaments, sheets and 
voids respectively. From left to right the columns in Fig. [6] 
show the distribution of angles formed between the principle 
axes, as mentioned above. 

The strength (or weakness) of the alignment signal can 
be quantified by <r, the average offset between a given distri- 
bution and a uniform one, calculated in terms of the Pois- 
son error. In practice, the average difference between the 
number of haloes found in a given bin and the number ex- 
pected from a uniform distribution is calculated in terms of 
the Poissonion error of a radnom distribution. If this is less 
than unity, then (on average) the distribution lies within the 
Poisson error of a uniform distribution. High a indicates a 
strong deviation from uniformity while a weak signal corre- 
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Figure 8. Same as Fig. [6] but for the angle formed between the 
halo's angular momentum vector J and the eigenvectors of the 
shear tensor. 



sponds to low a. In each panel we indicate the strength (or 
weakness) of the alignment signal by a how far away it is 
from a uniform distribution. Generally speaking we say the 
vectors are aligned if the signal is stronger than la. 

We begin by examining the full halo population ir- 
respective of web classification. There are are two salient 
points to be inferred from the top row of Fig. [6] 

(i) Halos tend to point parallel to £3 . While the distribu- 
tion of the angle between £2 and b is consistent with random 
both cos(f> = &3.a and cos0 = e\.c show statistically signif- 
icant parallel alignments. This implies that long axis of a 
halo tends to align itself with (=3. 
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Figure 9. The median of the fractional anisotropy of the velocity 
shear in the local environment of a halo as a function of halo 
mass. We present the mass dependence of FA for haloes found 
in knots, filaments, sheets and voids in black, green, blue and 
red respectively. Error bars represent the lcr spread about the 
median. Note that there exists a critical mass above which FA 
decreases with increasing halo mass. 



(ii) The strength of the halo alignment is mass dependent. 
Although the lowest mass haloes tend to show an alignment 
at just the ~20% level (for the angle between £3 and a, and 
ei and c), the more massive halo bins show considerably 
stronger alignments. This may reflect the tendency of more 
massive haloes to become more aspherical and for the weaker 
low mass alignment to be due to the ill-defined principle axes 
of low mass haloes (i.e. Allgood et al.|2 006). 



Turning to the second, third, fourth and bottom rows of 
Fig. [6] allows us to examine the web dependence of the shape 
alignment. Here too a mass dependence of the strength of 
the shape-alignment is seen for each web type. Surprisingly 
the alignment exists for all web types, implying that haloes 
are aligned with the shear tensor regardless of the web type 
they inhabit. That is, the main tracer of the halo alignment 
appears to be the principle axes of the shearing material, 
and not whether a region is expanding or collapsing along 
any of the principle axes. 

That said, the alignment for the low mass haloes ap- 
pears to slightly increase in filaments and sheets with re- 
spect to knots. Presumably this is due to the increased en- 
vironmental density in knots, which may have the effect of 
torquing low mass haloes and thus randomizing their align- 
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Figure 10. Same as Fig. [6] but for the angle formed between a 
subhalo's orbital angular momentum vector J and the eigenvec- 
tors of the shear tensor. 



merits. The alignment of intermediate and high mass haloes 
on the other hand shows a strengthening in voids (although 
low number statistics may play a role). 

We quantify the effect of web environment on the halo 
shape alignment (with I3) in the top panel of Fig. [7] where 
the median of the angular distribution as a function of web 
classification is plotted for the three halo mass bins. A uni- 
form distribution has a median of < |cos</>j >= 0.5 (shown as 
the dotted line). The visual interpretation of the alignment 
distributions is confirmed - the median angle in all mass bins 
varies by at most ~ 5° degrees for all mass bins in knots, 
filaments and sheets, and increases by ~ 15° degrees for the 



highest mass haloes in voids. The slight strengthening of 
the median alignment for low mass haloes in filaments and 
sheets (by around ~ 5°), is also visible. 

On the right most panel in Fig. [I] we show graphically 
the orientation of halo shape with respect to the large scale 
structure in a thin slice through the simulation. Each halo is 
represented by an ellipse in the x, y-plane whose axes ratio 
correspond to those calculated by the inertia tensor. While 
the statistical nature of the alignment signal is visible, it is 
strong enough to pick out by eye that haloes point along the 
structure they are located in. 



3.2.2 Alignments with respect to halo spin 

We now turn to the (kinematical) alignment of halo spin 
with the shear tensor. Here too, the absolute value of the 
angle formed between J and e* is presented since e% does 
not denote a proper geometric direction. 

In Fig. [8] we show the probability distribution of the 
angle formed between the angular momentum vector and ez, 
e.2 and e\ for all web types as well as for each individually. 
Let us begin be highlighting the similarities as found with 
the shape alignment in Sec. |3.2| and Fig. [5] 

The alignment between spin and the cosmic web is sig- 
nificantly weaker than that found using the shape in the pre- 
vious section. As in the shape-alignment, the parallel align- 
ment of spin with e.2 and perpendicular alignment with ei is 
web independent. Each web type (noisiness due to different 
sample size not withstanding) tells more or less the same 
story: the orientation and strength of halo spin alignment is 
mass dependent. 

The alignment between J and ez, is fairly well nuanced. 



As shown in previous work (e.g. Codis et al. 112012 Aragon 



Calvo et al. 2007 Trowland et al. 2012 I high mass haloes 
display perpendicular alignments with e.3., while lower mass 
haloes show a parallel alignment. Note however that in pre- 
vious work ez was defined specifically in the case of filaments 
or walls - here we show that the flipping of high mass halo 
alignment occurs with respect to ez regardless of web type. 



et al. 



2012 



The mass a t which the transition occurs (defined by Codis 
as 5 x 1O 12 M0) appears to decrease with web 
that halos in the intermediate mass bin show 
voids. This finding suggests that 



type, such 

a spin-alignment flip in 
spin alignments are qualitatively web independent but the 
strengths of these alignments are web dependent. 

The alignment of the halo spin axis with the other two 
eigenvectors (ei and §2) behaves more similarly with respect 
to mass and web type. In general, that spins tend to align 



with «=2 was first suggested by Navarro et al. (20041. Again, 
this alignment is strengthened for the most massive halo 
bin. Halo spin axes also tend to be perpendicular to ei, the 
direction of the greatest collapse (in knots, filaments and 
sheets) and of the weakest expansion in voids. 

As in the previous sub section, we quantify the strength 
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M < 10 115 10 11 ' 5 < M < 10 12 ' 5 M > 10 125 all Masses 



Total 8,742,899 132,658 17,771 8,893,328 
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8% 


3% 


17% 


8% 


Filaments 


23% 


27% 


57% 


23% 


Sheets 


35% 


47% 


25% 


35% 


Voids 


34% 


23% 


1% 


34% 



Table 1. The number of haloes in each mass bin and the fraction thereof in each web type, used in this study. Masses are given in units 

of r'M Q . 



of the alignment by looking at the median of each distribu- 
tion in the bottom panel of Fig. [7] A similar picture to that 
presented with respect to the long axis is seen, namely that 
the median alignment is roughly constant (changing by less 
than ~ 5°) for knots, filaments and sheets, and strengthens 
in voids for the two highest mass bins. That the high mass 
bin always has a median less than 0.5 reflects the tendency 
for the most massive haloes to orient themselves away from 
e-z- 



3.3 The effect of shear anisotropy on halo mass 

In the preceding section it was shown that alignments be- 
tween halo shape and the cosmic web appear dependent 
on mass and independent on web type (with the excep- 
tion of voids) . Although the halo spin- web alignment has a 
transitional mass above which the alignment tends towards 
perpendicular, this mass appears web dependent decreasing 
with decreasing web type (that is, higher in knots and lower 
in voids). As was already shown in Fig.jlJ each web type has 
a fairly different distribution of web sphericity, with voids 
defining the most aspherical structures while knots collapse 
more uniformly. Given that the transitional mass in the halo- 
spin alignment appears web dependent, it is interesting to 
ask if the FA, as measured by equation |6j of the velocity 
shear plays a role in the determination of halo mass. Or 
to put it another way: if the most massive haloes display 
the strongest alignment, is their environment (as defined by 
shear FA) responsible? 

In Fig. [9] we present the median shear FA as a function 
of halo mass for each halo subdivided into its web type. A 
number of interesting observations can be made. First, low 
mass haloes show no (or only a very weak) dependence on 
shear FA. At a specific halo mass however, the flat relation- 
ship disappears and an inverse relation between halo mass 
and shear FA is found. The most massive haloes exist in 
regions of near uniform collapse. Insight into halo assembly 
may be thus gained by this relation: a halo may only grow 
to large masses if the region it is in is uniformly collapsing. 
Large shears inhibit halo growth. On the other hand, smaller 
haloes can be born in a myriad of shear environments. 

This behavior is qualitatively web independent. How- 



ever it differs quantitatively when examining the relation 
for the four different web types. In a given web type, the 
mass at which the FA begins to become mass dependent 
differs - lower for voids, greater for filaments, sheets and 
knots. Although it is difficult to quantify the mass at which 
FA becomes mass dependent, if one examines Fig. [9] by eye, 
it appears to be similar to the transitional mass for the flip- 
ping of spin-alignment. It is thus implied that the FA of the 
halo's shearing environment determines halo mass and spin 
in a statistical sense. 



3.4 Subhalo alignments with respect to the 
cosmic web 

Since the web finder employed in this work probes the non- 
linear regions of the large scale structure, it is interesting to 
examine whether the positions or orbits of subhaloes deep 
within their parent potential are affected by the large scale 
structure. To do so, it is insightful to examine similar plots 
as shown in Fig. [6] and Fig. [8] but for subhalos, not host 
haloes. Although examining the spatial distribution of sub- 
haloes with respect to the LSS would be indeed very reveal- 
ing (e.g.[Libeskind et al.|2005||2007[|2009||Metz et al.|2008 



2009||Kroupa et al.||2005| |Wang et al.||2012| >, by looking at 



the ensemble of subhalo orbital angular momentum, spatial 
information is implicitely includeorl 

In Fig. [To] we present the distribution of the cosine of 
the angle formed between the orbital angular momentum of 
satellites and the three principle axes of the shear tensor. As 
in Fig. [6] and Fig. |8j we present the full distribution in the 
top row, and the alignment signal in knots, filaments, sheets 
and voids in the second, third, fourth and fifth row. 

The subhalo orbital alignment behaves differently than 
the halo spin alignment. Note that low mass host haloes will 
not contribute any subhaloes to this sample while massive 
host haloes will contribute many. We begin by examining 
the full signal irrespective of web type. 



3 This is simply because J sa t is, by definition, perpendicular to 
Tsat, the satellites position vector, and thus the distribution of 



Jsat reflects the distribution of r sa t- 
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Subhalo orbital angular momentum, denoted as J sat = 
fsat x u sa t, show a weaker (at the ~ 10% level) but still 
statistically significant parallel alignment with ei and per- 
pendicular alignment with £3. Since subhalos tend to orbit 
along (or in the plane of) their parent's major axis, and since 
the parent major axis tends to be aligned with S3, it follows 
that the sub halo orbits should tend to be perpendicular to 
£.3. Note that the flipping of subhalo orbital spin alignment 
is not seen for subhaloes when binned by parent mass. 

The strength of the alignment signal and its sense is 
quite clearly web and mass dependent. Fig. [lO] shows that 
the signal for a given parent halo mass bin gets weaker with 
web type from void to knot. Indeed the greatest mass bin 
shows a correlation with for all web types: ~ 7% in knots 
growing to ~ 15% in voids). For the intermediate mass bin a 
correlation only emerges in filaments, sheets and voids while 
in the lowest mass bin it is only seen in sheets and voids. 
The weakening of subhalo orbital spin alignment in knots is 
likely due to the denser environment where the abundance 
of neighboring haloes provides torques capable of weakening 
the correlation. Finally we note that in knots, the subhalo 
angular momentum alignment with ei is actually opposite 
in its sense to that in the other web types. 

In summary, the web environment and mass of each host 
halo has a direct impact on the alignment of its subhaloes. 
While only the most massive parents display an alignment 
of j S at with e3 in knots, even the least massive parent show 
the same alignment in voids. A similar picture is seen with 
respect to ei and £2. 



4 SUMMARY AND DISCUSSION 

The shape and spin of DM haloes and the angular momen- 
tum of their subhaloes define characteristic directions. These 
orientations are studied within the context of the large scale 
structure (LSS). The LSS exhibits a web-like structure, the 
so-called cosmic web, composed of knots, filaments, sheets 
and voids. These are characterized here by the velocity shear 
of the cosmic fluid. It is this velocity shear, engendered by 
the primordial perturbation field, which is believed to have 
endowed galaxies with their angular momentum through the 
process of tidal torquing. The main focus of this work is to 
study the alignment of the characteristic directions of each 
halo with those of the LSS, in order to attempt to under- 
stand the alignment of haloes with the LSS. 

The eigenvectors of the velocity shear tensor provide a 
common framework within which the directional properties 
of DM halos and of the LSS can be studied and compared 
( Hoffman et al.|2012||Libeskind et al.|2012[ ). The results pre- 
sented here are based on the study of the Bolshoi simulation 
(Klypin et al. 2011). The wide dynamical range and reso- 
lution of the simulation has enabled us to penetrate deeper 



into the non-linear regions surrounding haloes at the present 
time. 

DM halos and the cosmic web constitute two differ- 
ent manifestations of the dynamical evolution of the LSS. 
Both are scale dependent and have their specific directional 
properties. Any attempt to compare these two should start 
with selecting an effective dynamical resolution for both ha- 
los and the cosmic web. While a very coarse grid runs the 
risk of washing out any alignment signal (due to a Gaussian 
smoothing of the velocity field that is applied), a fine grid 
will probe the inner depths of each halo and not the large 
scale structure. Thus our selection of a grid whose cell size 
if ~lMpc//i, puts an upper limit of ~ 2 x W 14 M Q /h on the 
halo size we may probe. 

The eigenvectors of the velocity shear tensor provide an 
orthonormal basis of unit vectors that define the principal 
directions of the cosmic web. We use this basis to investi- 
gate the alignment of halo spin, shape (as defined by its 
moment of inertia) and subhalo orbital angular momentum 
with the respect to the LSS. While much previous work has 
been dedicated to examining these alignments, nearly all the 
literature focuses on just one or two web types individually 
since most cosmic web finders are tailor made for one spe- 
cific web type (e.g. filament or void finders). To the best of 
our knowledge, we are thus the first paper to examine halo 
alignments with respect to the velocity shear tensor, and in 
particular to all cosmic web elements simultaneously. 

Using this democratic approach, we have found a strong 
alignment of halo shape with the cosmic web. This is a 



confirmation of number of studies (e.g. |Altay et al. 2006 



Aragon-Calvo et al.|2007||Brunino et al.|2007| |Zhang ; et ah] 



2009) which looked at shape alignment of haloes in one or 



two individual web types. The alignment signal found here is 
statistically significant across the entire cosmic web and can 
be seen graphically in the right most panel on Fig. [I] This 
forms the bulk of one of our main results: the alignment be- 
tween halo shape and the cosmic web is web independent: the 
shear tensor traces this alignment. The alignment is visible 
for all web types when using the velocity shear to quan- 
tify the cosmic web. From this correlation, a conjecture may 
be made: the velocity shear is the sole determining factor 
of halo alignment. Further work at higher redshift is neces- 
sary to verify this conjecture. We defer this calculation to a 
future paper. Yet we still may generalize previous explana- 
tions which are limited by web type (i.e. accretion along a 
filament, etc) and suggest that the velocity shear is the more 
fundamental tracer of web independent halo alignments. 

The alignment of halo spin with the cosmic web has 
thus been examined in the context of previous findings which 
indicated that lower mass haloes tend to spin with their ro- 
tation axis aligned parallel to the large scale structure in 
filaments and walls. The alignment of high mass haloes in 
filaments on the other hand was found to tend to be perpen- 
dicular to the filament's axis. The mass where this flip occurs 
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was found in previous studies to be ~ 5 x 10 Mq ( Aragon 



Calvo et al. 20071. Studies such as Codis et al. (20121 at 



tributed this flip in the orientation of the rotation axis to the 
merger history - whereas low mass haloes grow by isotropic 
accretion of diffuse dark matter, higher mass haloes grow by 



the merger of subhaloes (e.g. Forero- Romero 2009 1 whose in- 
fall is correlated with the large scale structure ( Knebe et al. 



2004| |Zentner et al.|2005| |Libeskind et al.|2005||2011[ ). |Lac" 

erna & Padilla ( 2012 1 have suggested that the orientation of 
spins may be due to the assembly bias of high mass halos - 
a high mass halo dominates the accretion and hence the an- 
gular momentum in its local neighborhood at the expense of 
low mass haloes. Most of these studies implicitly assign the 
mass scale at which halo growth becomes merger dominated 
as the transitional spin flip mass of ~ 5 x 10 12 M, 



Trowland 



et al. (20121 have argued against this interpretation as they 



find that the spin axes of all haloes are oriented perpendicu- 
lar to filaments at high z, and it is just the spin axes of low 
mass haloes that evolve towards a parallel alignment. 

While considerably weaker than the shape alignment, 
we confirm the results of other studies that focus on the 
spin alignment. Yet we find that the mass at which the spin 
flips to be web dependent, hinting that instead of reflecting 
a physical process associated with the halo it is more likely 
a reflection of the interaction with the web environment and 
thus the scale against which it is measured. We have found 
that (at a fixed scale of ~lMpc//i) the transitional mass 
discriminating between perpendicular or parallel alignment 
decreases from ~ 10 12 ' 5 /t _1 M Q in knots to ~ 10 115 /i _1 M o 
in voids. Furthermore the mass at which the spin flips with 
respect to the large scale structure appears to be intimately 
related to the mass at which haloes become dependent on 
the shear's uniformity (see below). 

Since haloes are assigned the eigenvalues of the shear- 
ing material they are embedded in, we have compared this 
distribution with the full distribution of all cell eigenval- 
ues. These distributions differ indicating a significant bias 
for studies that attempt to reconstruct the large scale shear 



based on magnitude limited samples (e.g. Lee & Erdogdu 
20071. Just as galaxies represent a biased sample of the den- 
sity field (they trace high density peaks) , so too is the shear 
as measured by haloes a biased representation. 

The shear uniformity has been characterized using 
the fractional anisotropy FA (Basser 1995 see footnote 2) 
which measures how much of the shear magnitude is in an 
anisotropic element. A strong asymmetry in the distribu- 
tion is found, in the sense that knots appear to collapse more 
isotropically than voids expand (the collapse of a knot is not 
the expansion of a void played in reverse). Halo mass appears 
to correlate with FA above a certain mass: more massive 
haloes inhabit regions where the shear is more spherical. 
The mass at which this dependence occurs is roughly the 
mass scale for which spin-alignments flip. 

We have confirmed that the alignment of subhalo orbits 



is also well correlated with the large scale structure, albeit 
unlike the haloes the correlation is web dependent. Possibly 
owing to the denser environment in knots, the correlation 
is stronger in voids and weaker in knots. The large scale 
structure as defined by the velocity shear thus has a direct 
influence on the orbits of subhaloes deep in their parent po- 
tentials. This finding has a direct bearing on the correlated 



rotation of the Milky Way's own satellite system ( Metz et al. 
[20091 . 

As future surveys provide more and more detailed pho- 
tometry of galaxies in the framework of cosmic web, we will 
undoubtedly be able to infer more and more regarding the 
importance of environment as defined dynamically on the 
process of dm halo assembly and galaxy formation. Observa- 
tions of the spatial distribution of external satellite galaxies 
as well as the orientation of the spins of extra-galactic discs 
and spheroids with respect to the cosmic web, will help us 
understand the importance of the web in determining how 
galaxies acquire their angular momentum. 
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